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Damped Lyman-ct systems (DLAs) are useful probes of star formation and galaxy forma- 
tion at high-redshift (hereafter high-z). We study the physical properties of DLAs and 
their relationship to Lyman break galaxies (LBGs) using cosmological hydrodynamic 
simulations based on the concordance A cold dark matter model. Fundamental statistics 
such as global neutral hydrogen (H i) mass density, H I column density distribution func- 
tion, DLA rate-of-incidence and mean halo mass of DLAs are reproduced reasonably 
well by the simulations, but with some deviations that need to be understood better 
in the future. We discuss the feedback effects by supernovae and galactic winds on the 
DLA distribution. We also compute the [C II] emission from neutral gas in high-z galax- 
ies, and make predictions for the future observations by the Atacama Large Millimeter 
Array (ALMA) and Space Infrared Telescope for Cosmology and Astrophysics (SPICA). 
Agreement and disagreement between simulations and observations are discussed, as well 
as the future directions of our DLA research. 

Keywords: quasar absorption system; galaxy formation; cosmology. 



1. Introduction 

Recent studies suggest that our universe can be described well by a cosmological 
model called the concordance A cold dark matter (CDM) model with a matter 
energy density Q m = p m /p C rit — 0.3 and a dark energy density = Ph/ Pent — 
0.7. 1-3 Based on this backbone of structure formation, we would like to build a 
self-consistent model of galaxy formation that can describe the formation of disk, 
spheroid, and black holes. In particular, we would like to understand how the gas 
turned into stars as a function of cosmic time and environment. 

Observations of DLAs provide very unique opportunities for that purpose. Be- 
cause DLAs are observed in absorption, they are in principle an unbiased sample 
of H I in the universe, and provide us with information that are complementary to 
those obtained by direct observations of stellar emission. For example, with DLAs 

a DLAs are historically defined 4 to be quasar absorption systems with H I column density A r ni > 
2 X 10 20 cm -2 . See Ref. 5 for an extensive review on DLAs. 
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we can measure the global H I mass density in the universe as a function of rcdshift. 6 
It has been estimated that DLAs contain a significant fraction of H I gas in the uni- 
verse at z ~ 3, with the amount roughly equal to a half of current stellar mass 
density. 7-10 Therefore DLAs are considered to be important reservoirs of neutral 
gas for star formation in high-z universe, and their study would reveal the physical 
properties of neutral gas before they turn into stars. The conversion history of neu- 
tral gas into stars has been an active area of research over the past ten years, 11-16 
and the connection between DLAs and star-forming galaxies such as LBGs 12 is 
becoming one of the focus of current DLA research. 

In recent years, observers have made significant progress in constraining the 
statistical properties of DLAs using large samples of high-z quasars discovered by 
surveys such as the Sloan Digital Sky Survey (SDSS). 17 For example, Ref. 9 have 
extracted ss 500 DLAs from the SDSS Data Release 3 and tightened the constraints 
on the column density distribution function and the rate of incidence. The total 
number of detected DLAs has reached over 1000 as of early 2007 (Prochaska, private 
communication). Increased accuracy in these observables set tight constraints on the 
models of galaxy formation. 

In turn, theoretical understanding of DLAs remain unsatisfactory. Currently 
there are mainly two competing views on the physical origin of DLAs. One idea is 
that DLAs originate from rapidly rotating, massive disk galaxies. 18 This model can 
explain the large velocity widths of low-ionization metal absorption lines. However, 
Ref. 19 used hydrodynamic simulations to argue that the asymmetric profiles of low- 
ionization absorption lines can also be explained by the rotation, random motions, 
infall, and merging of protogalactic gas clumps equally well. Somewhat intermediate 
idea is intersecting multiple disks 20 or tidal tails observed in merging systems. 

There are a number of DLA studies using so-called semianalytic models of galaxy 
formation, 21,22 but it would also be desirable to study DLAs using cosmological hy- 
drodynamic simulations that directly solve the gas dynamics within the framework 
of ACDM universe. Ref. 23 & 24 focused on the effects of star formation and super- 
nova (SN) feedback on DLAs that were largely neglected by the earlier numerical 
work, 25-27 and showed that the distribution of DLAs could be significantly affected 
by SN feedback. Ref. 28 studied the effects of radiative transfer by post-processing 
the results of cosmological hydrodynamic simulations, but their simulations did not 
include models for star formation and SN feedback. It is still difficult to treat the 
radiative transfer dynamically from all galaxies in a large volume of space, as well as 
star formation and feedback by SNe and active galactic nuclei in large-scale cosmo- 
logical hydrodynamic simulations. This is primarily due to the heavy computational 
load required for such a simulation with full treatment of all physical processes. At 
the moment, we need to limit our calculation to a small number of objects in order 
to treat radiative transfer self-consistcntly with the star formation in simulations. 

In this review article, we summarize our recent numerical studies on DLAs using 
cosmological hydrodynamic simulations. We discuss the conflicts and agreement be- 
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tween simulations and observations, and the future directions of theoretical research 
on DLAs. 

2. Simulating DLAs 

We use cosmological smoothed particle hydrodynamics (SPH) code GADGET-2 
(Springel 2005) for our study. The code solves the hydrodynamics in a Lagrangian 
fashion by allowing the gas particles to cluster in high density regions, thereby 
providing higher resolution in galaxy-forming regions. The gravity is solved by a 
Tree-particle-mesh method that uses a tree algorithm for the short-range force 
and particle-mesh technique for the long-range force. It also adopts the entropy- 
conservative formulation 29 to remedy the overcooling problem, which previous gen- 
eration of cosmological SPH codes suffered from. The code includes radiative cool- 
ing, heating by a uniform UV background radiation of a modified Haardt & Madau 
spectrum, 30-32 star formation, supernova feedback, galactic wind, 33 and a sub- 
resolution multiphase model of interstellar medium (ISM). 34 

We employ a series of simulations with different resolution, box sizes (from 
Lbox = 4/i _1 Mpc to 100/i -1 Mpc) and feedback strengths to study their effects. 
One of the focus of our study is to examine the effects of galactic wind feedback 
on the distribution of DLAs. In this article, we summarize our results by refer- 
ring to two cases: 'no wind' feedback and 'strong wind' feedback, corresponding 
to the wind speed of v w = and 484 km s -1 . See Ref. 24 for more details of our 
simulations. All of our simulations are based on a concordance ACDM cosmology 
with cosmological parameters (O m , ^b, 08, h) = (0.3,0.7,0.04,0.9,0.7), where 
h = ffoAlOOkms^Mpc- 1 ). 

3. Fundamental DLA Statistics 

3.1. HI Column Density Distribution 

The Hi column density distribution function /(iVgi) is the most fundamental statis- 
tic for any Hi quasar absorption systems. It counts the number of columns in bins 
of [Am, Am + dNm] per unit absorption distance dX, i.e., dN — f (Nm)dNBidX , 
where dX = ^(1 + z) 3 cdt — (1 + z) 2 dr. The quantity dN can also be inter- 
preted as an area covering fraction in the sky along the line element cdt. 

Ref. 23 examined /(Ahi) in a series of cosmological SPH simulations with dif- 
ferent box sizes and resolution. The left panel of Fig. [T] shows the results for runs 
with no wind and strong wind feedback, together with the observed data points. 9 At 
high column densities (log Am > 21), the agreement between the strong feedback 
run and observed data is quite good. The run with no wind feedback significantly 
overpredicts /(Ahi)- This is because the galactic wind feedback reduces the num- 
ber of columns with large Am by ejecting the gas from star-forming regions and 
heat the Hi gas. This agreement at log Am > 21 for the strong feedback run is 
encouraging and shows the effectiveness of our current wind feedback model. One 
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Fig. 1. Le/t: Column density distribution function /(-/Vhi) for the runs with no wind (open 
squares) and strong wind feedback (solid triangles). The red short-dashed line is the gamma 
function fit to the observed data points. 9 Right: Comoving DLA cross section as a function of 
halo mass. The shaded contours show the distribution of halos for the run with strong wind feed- 
back. The open triangle symbols are the median points in each mass bin, and the solid line is the 
least-square fit to the median points. The dashed line is the same fit for the run with no wind 
feedback. 



might speculate that the H I in the overpredicted columns in the no feedback run 
could turn into molecular hydrogen, therefore the apparent disagreement between 
the simulation and the data points is not a serious problem. However, Ref. 35 showed 
that the H2 column density distribution forms a natural extension to f(Nm) for 
local galaxies, and that the connection point between the two distribution is at 
(log ./Vhi, log /(iVm)) ~ (22, —26). Since the degree of overprediction in the no feed- 
back run is more than a factor of 2 at logiVni < 22, it seems difficult to hide all 
of the overpredicted H 1 in the form of H2 if similar atomic and molecular physics 
hold in high- 2: galaxies as in local galaxies. 

Our simulations underpredict the number of columns with 20 < log TVjji < 21 
regardless of the strength of wind feedback. Physical reasons for this underestimate 
is not clear at this point. Several test runs show that this feature does not change 
very much with varying threshold density or time-scale for star formation, or varying 
wind feedback strength. Any combination of the following could be the cause for 
this underestimate: inadequate resolution, angular momentum transfer problem, 
inadequate cooling due to lack of metal cooling in the current code, too simplistic 
models for star formation and SN feedback, or poor treatment of radiative transfer. 
We plan to address these issues in our future work. 
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3.2. Neutral Gas Mass Density 

For a given observed path length cAt, the H I mass density within the survey volume 
can be estimated as follows (in proper units): 

Phi = ^ / NmaX N m f(N m )dN m AX = g ° WH(1 + z)3 N m f(N m )dN m , 

(1) 

where mn is the mass of a hydrogen atom, and [-/V m ; n , N max ] is the observed range 
of iVni- If the observed range of Ahi is broad enough and covers up to sufficiently 
high iVm values, then one can expect that the calculated quantity is representative 
of the true value in the universe, and the mass density parameter of neutral gas can 
be estimated as r2 n cutrai = MPHi/Pcrit, where \i « 1.3 is inserted to take helium into 
account. It is known observationally 7-9 that DLAs dominate the H I mass density at 
z ~ 3, with the amount roughly equal to half of the present-day stellar mass density. 
This suggests that DLAs are important reservoirs of neutral gas for star formation 
in high-z galaxies, and that the neutral gas in DLAs at z — 3 are converted into stars 
between z = 3 and z = 0. Ref. 23 showed that the observed estimates of f2 nC utrai 
can be bracketed by the simulation with no wind and strong wind feedback. It can 
also be shown that high-A^ji systems dominate flm, therefore it is important to 
include a fair number of high-TVHi systems in the sample to obtain an accurate fit 
to f(Nm) and estimate f^m reliably. 



3.3. DLA Cross Section vs. Halo Mass 

It is useful to quantify the cross section of DLAs (ctdla) as a function of dark 
matter halo mass, because we have a theoretical framework to compute the halo 
mass function. 36,37 Observers can identify DLAs only in the quasar sight- lines, but 
there could be multiple neutral gas clouds within the same halo that do not cross 
the quasar sight-line. Therefore it is theoretically more plausible to characterize the 
occurrence of DLAs as a cross section rather than the number of clouds in each 
halo. 

Ref. 23 estimated the relationship between the total DLA cross section crg 5 LA (in 
units of comoving h~ 2 kpc 2 ) and the dark matter halo mass (in units of /i _1 M Q ) at 
z = 3 as 

log ^dla = <* (log M halo - 12) + p, (2) 

with slopes a = 0.7—1.0 and normalization (i = 3.9 — 4.2, depending on the 
numerical resolution and the strength of galactic wind feedback. The slope a is 
always positive and massive halos have larger DLA cross sections, but they are 
more scarce compared to less massive halos. Two important features are: (1) As the 
strength of galactic wind feedback increases, the slope a become steeper while the 
normalization /? remains roughly constant. This is because a stronger wind reduces 
the amount of neutral gas in low-mass halos at a higher rate by ejecting the gas out 
of the potential well of the halo. (2) As the numerical resolution is improved, both 
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the slope and the normalization increase. This is because with higher resolution, 
simulations can resolve higher densities, and star formation in low-mass halos can 
be described better. As a result the neutral gas content is decreased due to winds. 
On the other hand, a lower resolution run misses the early generation of halos and 
the neutral gas in them, resulting in lower ^hi at high-z. 

Earlier numerical work 25-27 largely neglected the effects of star formation and 
wind feedback, therefore they often overestimated the DLA cross sections in low- 
mass halos, resulting in higher value of predicted rate-of-incidence (see § 13. 4| . In 
order to match their results with observed data, they introduced a minimum halo 
mass that DLAs populate. Our simulations show that star formation and feedback 
affect the DLA distribution, therefore, we need to study their effects further with 
better resolution and improved models of star formation and feedback. 



dr r°° 

-(> M) = — / n dm (M') < LA (M') dM', (3) 
az J M 



3.4. DLA Rate-of -Incidence 

The probability of finding a DLA for a given observed path-length can be ex- 
pressed as a rate-of-incidence. Theoretically, we can compute the cumulative rate- 
of-incidence of DLAs per unit redshift above certain minimum halo mass M as 

dJV D LA , 

, M) = — 

where n<i m (M) is the dark matter halo mass function and dr/dz = c/H(z) with 
H(z) = i/o\/fim(l + z ) 3 + Cl\ for a flat universe. In order to carry out this integral, 
Eq. © can be used. 

The left panel of Fig. [2] shows that the run with no wind feedback agrees with 
the observed rate if we sample the halos down to Mhaio — 1O 8 5 /i _1 M . The run 
with strong feedback underestimates the total DLA rate-of-incidence by about 0.3 
dex due to the underestimate of f(Nm) at logA^i < 21. Based on the results of 
hydro simulations, Ref. 38 derived the mean scaling relation between halo mass and 
galaxy i?-band AB magnitude as 

Rab = - 2 • 5 log M ha i + C - 5 log h 7a , (4) 

where C = 55.03 (no feedback) and 57.03 (strong feedback). This relation was used 
to convert the halo mass into apparent Rab magnitude in the left panel of Fig. O 
The arrows in the figure indicate that the wind feedback reduces the contribution 
from low-mass halos to the DLA rate-of-incidence. 

One can also calculate the differential probability distribution of DLAs as a 
function of halo mass as 

dAl "- A =^[Mn(M)Hl0)}a^ LA (M). (5) 



dz d log M dz 

See Ref. 38 for the derivation of Eq. and ([5]) . Since stronger wind feedback sup- 
presses the contribution from low-mass halos to the total DLA rate-of-incidence, the 
distribution shifts toward higher mass halos as the feedback strength is increased, 
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Fig. 2. Left: Cumulative DLA rate-of-incidence for all galaxies above certain magnitude limit. An 
observational estimate 9 is shown by the hatched band. Right: Probability distribution function of 
DLA rate-of-incidence as a function of dark matter halo mass. As the feedback strength increases, 
contribution from low-mass halos decreases, and the distribution shifts toward higher masses with 
a peak emerging at M^aio — 10 12 h — 1 Mq. 

and eventually a peak in the distribution appears at Mhaio — 1O 12 /i _1 M0 as shown 
in the right panel of Fig. O This mass-scale is interesting, because recent observa- 
tional studies 39 ' 40 of cross-correlation between DLAs and LBGs suggested similar 
characteristic halo masses for DLAs. Ref. 41 also estimated the mean halo mass 
of DLAs using numerical simulations (see Ref. 38 for a detailed comparison of our 
results with theirs.). 

3.5. Impact Parameter Distribution 

Another important way to characterize the spatial distribution of DLAs is to exam- 
ine the distribution of impact parameter "6" between DLAs and the nearest galaxy. 
For a given DLA sight-line in the simulation, we search for the nearest galaxy on 
the projected plane, and determine the ^-parameter for each DLA. As the left panel 
of Fig. [3] shows, the ^-parameter distribution is quite narrow in the current simu- 
lations, and over 80% of DLAs have a galaxy within physical kpc. This large 
fraction of DLAs with small ^-parameter is due to the large number of low-mass 
halos in a CDM universe and DLAs associated with them. As the Figs. 1 — 3 in 
Ref. 24 show, DLAs in lower mass halos are very compact with physical sizes of 
a few kpc, and coincide with star-forming regions very well. One might consider 
that this is plausible given the fact that the physical radii of LBGs are 2 — 3 kpc, 42 
however, DLAs could also originate from extended Hi clouds observed as, e.g., tidal 
tails of merging systems. We do not seem to have such extended DLAs associated 
with low-mass halos in our current simulations. 
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Fig. 3. Left: Normalized impact parameter distribution for runs with and without galactic wind 
feedback. The results with limited search for the nearest galaxy to Rab < 28 mag are also shown. 
Right: Cumulative comoving number density of DLAs as a function of apparent Rab magnitude 
for three different values of assumed DLA physical radius rrjLA = 1^70 kpc, S/i^kpc, and 



20 h 7 } kpc. The covering area of each DLA is Afjla 



For each value of rrjLAi results 



from no feedback run (long-dashed line) and strong feedback run (solid line) are shown. The thick 
dot-dashed line is the cumulative comoving number density of LBGs obtained by integrating the 
observed luminosity function. 46 The data point at Rab = 25.5 shows the observed comoving 
number density of LBGs JVlbg = 4 x 10~ 3 h~ 3 Mpc -3 . 



If we limit the search for nearby galaxies to those brighter than Rab — 28 mag, 
then the distribution becomes broader and roughly 80% of DLAs have a galaxy 
within physical lOh^Q kpc. If we increase the strength of galactic wind feedback, 
then we reduce the fraction of DLAs in low-mass halos, and a larger fraction of 
DLAs will be in more massive halos that are farther away, resulting in a broader 
6-parameter distribution. It is not clear at this point whether these conclusions 
from the current simulations are consistent with current observations given the 
limited sample. 43,44 When a larger sample of DLA host galaxies is constructed, the 
comparison will become more interesting. 

We can also constrain the spatial distribution of DLAs using the cross-correlation 
function 45 between DLAs and galaxies. Ref. 39 computed the cross correlation func- 
tion using 211 LBGs and 11 DLAs, and showed that there is a strong correlation 
between the two population, comparable to the LBG auto-correlation. This suggests 
that many DLAs are physically related to LBGs and quite possibly in the same dark 
matter halos as LBGs. We will present the DLA-LBG cross-correlation function in 
our simulations in the future. 
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3.6. Number Density of DLAs 

It would be useful to quantify the number density of DLAs similarly to the galaxy 
number density. However doing so is not easy both observationally and in numerical 
simulations, because DLA gas clouds are not detected in emission and they are often 
extended, connected, and exhibit very complicated morphologies as shown in Fig. 1 



Theoretically, it is possible to calculate the mean number density of DLA clouds 
as follows. Assuming that the average physical covering area of each DLA gas cloud 
is fixed with Aula — "" r DLA' then the cumulative number density of DLAs above 
a certain halo mass M is 



where nam{M) is the dark matter halo mass function, and a^^(M) is the physical 
DLA cross section as a function of halo mass. Note that ctdla is the total DLA cross 
section of each dark matter halo; in other words, if there are 100 DLA clouds in 
a massive halo, then this halo has a total DLA cross section of cr^j^ = IOOAdla- 
Then, using Eq. (fj[|) and © , we calculate the cumulative comoving number density of 
DLAs as a function of apparent Rab magnitude Adla(< Rab) as shown in the right 
panel of Fig. [3] This figure shows that the number density of DLA clouds is higher 
than that of LBGs if their physical radius is smaller than tola = 20h^Q kpc. Based 
on a similar calculation, Ref. 47 suggested possible relationship between outflows 
from LBGs and DLAs, although this idea is yet to be confirmed observationally. 

3.7. Metallicity, TVhi, and Star Formation Rate 

Measurements of gas metallicity in DLAs and galaxies provide useful information 
on star formation and chemical enrichment history. The mean metallicity of DLAs 
is known 48 to be w (1/30)Z Q , and that of LBGs to be « (1/3)Z . 

In Figure HI we show the distribution of gas metallicity as functions of A/hi and 
projected star formation rate Ssfr for the two strong wind feedback runs (Q5 and 
G5 run) at z — 3. The G5 run has a broader distribution than the Q5 run because 
of a larger box size (Lbox = 100/i _1 Mpc for the G5 and 10/i _1 Mpc for the Q5 
run), but lower median metallicity at lower TVhi values due to lower resolution. In 
the left panel, we see that there are no observed DLAs 49 with log Am > 22.0 and 
log(Z/Z ) > —1.0, while such sight-lines exist in our simulations. These columns 
with high Ahi and high metallicity are probably associated with central cores of 
LBGs that are actively forming stars. The occurrence of such columns are very 
rare as can be seen from the rapidly declining /(Am) function, and they have 
very small cross-sections. There are other ideas to explain the absence of DLAs 
at log Am > 22.0. One popular idea is the extinction of background quasars by 
dust. 50,51 Another idea is that the conversion of Hi to H2 molecule determines the 
maximum value of Am- 52 It is not fully clear at this point which of the above three 



of Ref. 24. 
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Fig. 4. Left: Distribution of DLA columns on the plane of Hi column density vs. metal- 
licity are shown for the Q5 run (L^ox = lO/i^Mpc; magenta contour) and the G5 run 
(Lbox = 100/i — 1 Mpc; black contour) at z = 3 with the median points for each logTVjji bin. 
The blue open triangles are observed data. 49 Right: Same as the left panel for the plane of pro- 
jected star formation rate vs. metallicity. The skewed contour shape at logEgpH, ~ —3 for the G5 
run is a resolution effect. The blue data points are observational estimates by the C II* technique, 54 
and the open square symbol represents LBGs roughly. 



explanations (or any combinations ol them) is the correct one. At log Am < 22.0, 
the Q5 run might be overpredicting the median metallicity of DLAs compared to 
observations. 

The right panel of Fig. [4] shows a similar distribution as the left panel. This is 
because Ssfr and Ahi are correlated with each other strongly (i.e., Ssfr oc A 1 ^ 4 , 
the Kennicutt Law; 53 see also Fig. 8 of Ref. 24), and metallicity also correlates with 
star formation positively. The metallicity for the columns with log £sfr ^ in the 
simulation might be slightly too high compared to the observational estimates of 
LBGs. The observational estimates by the C II* technique 54 (blue data points; only 
positive detections are shown) nicely fall onto the simulated contour. 

4. Predicting the [C II] Line Emission 

The [Cn] emission line at A = 157.74 /im originates from the 2 P3/2 P\ii fine 
structure transition of C + . This emission line is often the brightest emission in 
the spectrum of galaxies, and it is suggested to be the dominant coolant of diffuse 
ISM at temperatures T < 5000 K. 55-58 Therefore, it has been suggested that the 
[C II] emission is potentially another method for finding cold neutral gas in high-z 
universe. 59,60 

Locally, [Cn] emission has been detected from both 1) dense star- forming gas 
irradiated by UV radiation from young stars 61-64 and 2) extended, diffuse cold 
components of ISM in quiescent late-type galaxies. 65-68 At high-z, there is only a 
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few tentative detections 69, 70 of [Cn] emission from quasars at z = 6.42 and z = 4.7, 
but none for more normal (i.e., not quasars or active galactic nuclei) star-forming 
galaxies such as LBGs. Detections of [Cll] emission from LBGs, DRGs, 71 ' 72 and 
star-forming BzK 73 galaxies might become possible in the near future by ALMA 
and SPICA. 

We might be able to set important constraints on galaxy formation models using 
[C II] detections. In particular interferometric maps of [C n] emission would reveal 
the spatial dimensions of neutral gas that is responsible for DLAs, a property that 
has so far eluded detection. By combining these measurements with the velocity 
widths of the [C n] emission lines, it will be possible to infer the dark matter halo 
masses of the sources, and the [C n] luminosity would tell us the heating rate of the 
gas. 

Ref. 74 attempted to make a realistic estimate of [Cn] emission from high-;? 
galaxies by coupling the analytic model of Ref. 75 and cosmological hydrodynamic 
simulations. We first used the star formation rate and metallicity of gas in our 
simulations to estimate the amount of cold neutral medium (CNM; typically T ~ 
100 K and n ~ 10 cm -3 ) in each halo, and then computed the [Cn] luminosity by 
multiplying the [C n] luminosity per H I atom estimated by the analytic model to 
the CNM mass. As shown in the left panel of Fig. El we find that the [C n] luminosity 
can be described roughly by 



where C\ = 10 41 ergs -1 (no wind) and 10 40 ergs -1 (strong wind feedback). The 
strong wind feedback reduces the amount of CNM by ejecting the gas from star- 
forming regions, resulting in lower [C n] luminosity by a factor of ~ 10 compared 
to the no feedback run. 

The right panel of Fig. [5] shows the probability distribution function of [C il] 
sources at z — 3 as a function of [C n] flux density. For the runs with strong 
feedback, we see that only a few percent of all the sources would have [C n] flux 
brighter than 1 mJy Most of the [C n] sources will be very faint, low-mass galaxies 
with S v < 0.1 mJy due to the large number of low-mass halos in a ACDM model. 
The run with no wind feedback gives more optimistic result in terms of the fraction 
of bright sources, suggesting that ~ 30% of [Cn] sources have S v > 0.01 mJy at 
z = 3. These predictions can be tested by ALMA and SPICA in the future. Even 
though the [Cn] observation of normal high-z galaxies is currently difficult, the 
reward of [C n] detection from such sources would be quite significant, because such 
measurements could directly constrain the amount of CNM and the properties of 
ambient ISM, which is otherwise difficult to achieve. 

The crude calculation by Ref. 74 was satisfactory as a first step, but we need to 
treat the radiation from local sources and the ionization of gas more accurately. To 
this end, we plan to use the photoionization code CLOUDY 76 and refine our [C n] 
calculations in the future. 




(7) 
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Fig. 5. Left: [Co] luminosity as a function of halo mass at z = 3. The three sets of points are 
all from strong wind feedback runs: Q5, D5, and G6 runs with box sizes Lt, ox = 10, 34, and 
100/i — 1 Mpc, respectively. Right: Cumulative probability distribution function of [Cn] sources at 
z = 3 as a function of flux density. The result from the run with no wind feedback (the 03 run 
with L box = 10h~ 1 Mpc, long-dashed line) is also shown. The two magenta dot-dashed lines are 
obtained 74 by integrating the observed luminosity function 46 of LBGs. The three vertical dotted 
lines indicate the flux density limits of Sv = 0.1,0.01, and 0.001 mjy. the diagram. 



5. Conclusions & Future Research 

In this article, we reviewed our DLA research using cosmological SPH simulations. 
We achieved some successes in explaining the global DLA statistics crudely, but 
there are details that need to be studied further. In particular, our simulations 
underpredict /(iVgi) at logA^i < 21.0. To cure this problem, we need to make 
improvements on the following points: 1) increase mass and spatial resolution, 2) 
treat radiative transfer from local radiation sources and photoionization of gas more 
accurately, 3) improve the models for star formation and feedback. Also, our simu- 
lations are probably still suffering from the angular momentum transfer 77 problem, 
which we hope to remedy as we work on the above three points. 

The possible connection between DLAs and LBGs is quite interesting. Ref. 75 
detected the C II* A1335.7 absorption in about a half of randomly selected DLAs, 
and showed that local heat sources (such as LBGs) are required to match the heating 
rate inferred from the [Cn] cooling rate. Their work and Ref. 39 suggest close 
connection between DLAs and LBGs. It would be natural for the two population 
to be closely related to each other, considering the fact that DLAs are important 
reservoirs of neutral gas for star formation at high-z, and that LBGs contribute 
significantly to the total star formation rate density in high-z universe. 

Just like the multi-wavelength observations of galaxies can reveal different prop- 
erties of galaxies, we can also tackle the mysteries of DLAs using different techniques 
including H I absorption, C n* and other metal line absorption, and [C n] emission. 



I 



I 
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Through continuous collaboration between observation and theory, we hope to clar- 
ify the physical nature of DLAs in the future. 
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